Phase Separation close to the density-driven Mott transition in the Hubbard-Holstein 

model 
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The density driven Mott transition is studied by means of Dynamical Mean-Field Theory in the 
Hubbard-Holstein model, where the Hubbard term leading to the Mott transition is supplemented by 
an electron-phonon (e-ph) term. We show that an intermediate e-ph coupling leads to a first-order 
transition at T — 0, which is accompanied by phase separation between a metal and an insulator. 
The compressibility in the metallic phase is substantially enhanced. At quite larger values of the 
coupling a polaronic phase emerges coexisting with a non-polaronic metal. 
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Since long time it is known that strong electron- 
electron (e-e) interactions can drive a metallic system 
insulating [lj. The correlation caused Metal-Insulator 
transition (MIT), usually known as the Mott transition, 
occurs in several compounds like V2O3, CaTiOz and is 
often theoretically investigated in the framework of the 
Hubbard Hamiltonian. However, although this latter cer- 
tainly captures the fundamental properties of correlated 
electronic systems, it lacks the lattice degrees of freedom, 
which can play a crucial role in many respects. Therefore, 
not only the interplay between the electronic correlation 
and the electron-phonon (e-ph) interaction is a relevant 
issue in general 0, Q , but interesting specific effects can 
arise from it, like in doped fullerenes, where a substantial 
e-ph coupling associated to Jahn- Teller distortion can be 
synergetic with the large electron-electron interaction to 
produce the observed rather large superconducting tem- 
peratures p|. Many experimental results also point to- 
ward an important role of phonons in the strongly cor- 
related cuprates. In this light, it has also been proposed 
that the e-ph coupling could induce charge instabilities 
in the cuprates [2, la] thereby relating the interplay be- 
tween strong e-e and e-ph interactions to the issue of 
electronic phase separation (PS), stripe formation, and to 
the isotopic dependence of the pseudogap crossover tem- 
peratures in these materials . This connection between 
e-ph coupling and charge instabilities suggest that lattice 
degrees of freedom may greatly influence the thermody- 
namical stability of a strongly correlated electron system. 
The present work precisely focuses on the role of the e-ph 
coupling in modifying or emphasizing the physical prop- 
erties of electrons close to their Mott insulating status. 
Specifically we study the role of phonons in changing the 
order of the density-driven MIT and in inducing PS close 
to the transition. 

From the technical point of view we exploit the great 
progress made in this field with the dynamical mean-field 
theory (DMFT), which emerged as the first theoretical 
tool able to give a complete characterization of the Mott 



transition^- This technique, which becomes exact in the 
infinite coordination limit 8] , is in fact able to treat with 
the same accuracy both the metallic and the insulating 
phase, and it has been successfully applied to the Mott 
transition in the Hubbard model. Moreover, DMFT al- 
lows to treat on the same footing the e-e and e-ph inter- 
actions, which is crucial to our study. The main draw- 
back of this approach is the momentum-independence of 
the self-energy, and the neglect of non-local interactions. 
This is not expected to introduce significant biases when 
the interactions are local, like in the model we study. 

One of the most striking features of the DMFT treat- 
ment of the Mott transition is that both the interaction- 
driven transition at half- filling 0, [^} and the density- 
driven transition I ll| are characterized by a coexistence 
of solution in some region of parameters. In the latter 
case, which we address in this paper, there is a region of 
chemical potentials, delimited by two curves fi c i(U) and 
Hc2(U), in which a metallic solution with n 7^ 1 and an 
insulating one with n = 1 coexist Naturally, the sta- 
ble solution is the one that minimizes the grand-canonical 
free energy. At T = the metal is stable in the whole 
coexistence region, and is continuously connected to the 
insulator ^lj. Moreover, the metal joins onto the in- 
sulator with a finite slope in the n-/z cu rve, i.e., it has 
a finite compressibility n = dn/d^i [ill fl2j |. At finite 
temperature the energetic balance gets more involved, 
and a first-order transition occurs along an intermediate 
line between /z c i and /i C 2, leading to PS between the two 
phases 0. 

In this work we study how this scenario is modified by 
the inclusion of a short-range interaction with phonons. 
We consider the e-ph interaction in its simplest model re- 
alization, the Holstein molecular crystal model, in which 
tight-binding electrons interact with local modes of con- 
stant frequency [Til]. The Hamiltonian is 
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where Cj (cj) and a, (aj) are, respectively, destruction 
(creation) operators for fermions and for local vibrations 
of frequency u> on site i, t is the hopping amplitude, g 
is an e-ph coupling. 

Previous studies have identified two different dimen- 
sionless e-ph couplings, A = 2g 2 /u>ot, which measures 
the energetic convenience to form a bound state, and 
a = g/u>o which controls the number of excited phonons. 
The relevance of each coupling depends on the adiabatic 
ratio 7 = LUo/t: If 7 is small A is the control parame- 
ter, while for large 7 the physics is mainly controlled by 
a[ljj]. Here we work with 7 = 0.2, and therefore we use 
A to measure the strength of the e-ph coupling. 

In DMFT, the lattice model is mapped onto an im- 
purity problem subject to a self-consistency condition, 
which contains all the information about the lattice. In 
the case of the Hubbard-Holstein model we have an An- 
derson impurity with a local phonon on the impurity site. 
If we work in the Bcthe lattice of half-bandwidth t the 
self-consistency enforcing the DMFT solution is given by 

4 iuj n - e k 

where and Vu are the energies and the hybridization 
parameters of the Anderson impurity model. We use 
exact diagonalization to solve the impurity model 
The technique consists in restricting the sum in Eq. (J2J to 
a finite and small number of levels N s — 1 . The discretized 
model is then solved using the Lanczos method, which 
allows to compute the Green's function at T = 0. The 
method converges rapidly as a function of N s , and just 
a few levels are enough to achieve convergence. Here we 
mainly present results for N s — 10, having checked in 
specific cases that no measurable change occurs in the 
physical quantities by increasing AT.[l7|. 

As anticipated above, we discuss how the density- 
driven Mott transition is affected by the e-ph interaction. 
The main trend, which can be anticipated on intuitive 
grounds, is that the e-ph interaction favors the insulat- 
ing phase with respect to the metal. We note that in Ref. 
lid the same effect is obtained by raising the tempera- 
ture, because in that case the insulating state has a larger 
entropy. Here, when the e-ph interaction is turned on in 
the metallic state, two competing effects take place. On 
one side, the electron gains potential energy if the lat- 
tice is distorted, while, on the other side, the effective 
mass is enhanced, leading to a loss in kinetic energy. As 
a result of this competition, the electrons can not com- 
pletely exploit the e-ph coupling to lower the energy. On 
the other hand, the electrons in the Mott insulator are 
already localized, so that they can gain potential energy 
by coupling with phonons without losing kinetic energy. 
This simple argument already tells us that the energetic 
balance between the two phases will be changed in favor 
of the insulator. For the same reasons we can also expect 



the region of \i where the insulator exists to be larger, 
i.e., that fx c i decreases with increasing A [lij ]. 

We notice that simplified approaches, as the one based 
on Lang-Firsov and squeezing transformations ^IJ also 
suggest that the insulating behavior may be favored by 
the electron-phonon interaction. Within this approach, 
the Hubbard-Holstein model is transformed into a Hub- 
bard model with reduced hopping, in which the insulating 
behavior is clearly favored with respect to the case where 
the electron-phonon interaction is absent. 

These naive arguments are confirmed and based on 
solid ground by the DMFT calculations. We have studied 
in detail both the solutions as a function of A for values 
of U/t larger than the critical value U/t ~ 3 for the Mott 
transition at half-filling 0. In Fig □ we plot the density 
as a function of the chemical potential for U/t — 10 and 
various values of A. For A = (not shown) we recover the 
known results for the Hubbard model, with two solutions 
in a really small /i interval, and the metal is stable in 
the whole coexistence region, as it can be checked by 
comparing the grand canonical potentials = E — fj,N 
of the two phase. In the insets of the various panels we 
plot the difference O met — f2j ns . For small values of A 
we observe a slight increase of the coexistence region, 
but the metal stays at lower energy than the insulator 
(see upper panel of Fig. QJ. For A > 1 the two curves 
start to cross and the grand canonical potential of the 
insulator becomes lower than the one of the metal for 
some range of chemical potential (in the middle panel 
we show the result for A = 2), leading to a jump of the 
occupation number by continuously varying the chemical 
potential. The transition becomes therefore of first order, 
and the chemical potential has a plateaux with respect 
to the density. As a result, the system is not stable in 
the interval of density between which the jump occurs, 
and it undergoes PS between the insulating solution with 
n = 1 and a metallic solution with some density n ps (X). 
n ps is found to be a decreasing function of A. In Fig. [21 
we show the PS region for U/t = 10. 

If we further increase the coupling, the effect of the 
e-ph interaction becomes even more dramatic. In the 
case of U/t = 10, and A ~ 2.2, the insulating solution 
continuously evolves into a phase with density slightly 
smaller than 1. Also for this coupling the insulator is 
energetically favored in the coexistence region (see in- 
stet of lower panel in Fig. QJ. An inspection to other 
physical quantities allows us to better characterize this 
strong-coupling phase as a polaronic state. In fact, if we 
compute the phonon displacement distribution function 
P(X) = {X\ipo)( , ipo\X), the evolution from the n = 1 
insulator to this state is accompanied by a change from 
a unimodal distribution of the phonon displacement to a 
bimodal distribution characteristic of the polaronic state. 
On the other hand, the metallic solution does not display 
polaronic features for these values of A and so the PS in 
this region, takes place between a polaronic insulator and 
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FIG. 1: The n-^i curves for U/t = 10 for A = 0.6, 2, 2.2. Solid 
circles joined by a dashed line denote the metallic solutions, 
and the solid line marks the insulating one. The insets show 
the difference between the grand canonical potentials of the 
two phases Afl — Q me t — ^ms- In the two lower panels this 
quantity is negative signaling the stability of the insulator. 
An arrow marks the transition point on the n-jj, curves 

Our DMFT results show that the Hubbard-Holstein 
model displays PS close to half-filling, due to the first- 
order phase transition between the insulator and the 
metal. Even if PS occurs close to half-filling, just like in 
the large— N treatment of the same model, the present 
result is physically quite different. In Ref. [2(| a diver- 
gent compressibility was indeed found within the metallic 
phase. Here this instability is prevented by the stabiliza- 
tion of the insulator. Nevertheless, we can follow the 
metallic solution in the regime where it is metastable, 
and check whether the e-ph interaction may lead to a 
divergence, or at least to a sizeable enhancement of the 
charge compressibility. 

In the infinite-?/ limit, the large— TV approach has 
shown that the charge compressibility is well represented 
by a simple RPA-like formula, 
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where k(A) is the compressibility for a given value of 
the e-ph coupling A, TV* is the quasiparticle density of 
states per spin at the Fermi level, (Fq ) e (ph) is the elec- 
tronic (phononic) contribution to the symmetric Lan- 
dau amplitude, and (-F(f)p/, = —AN*g 2 /uo , which dif- 
fers from our A for the presence of the quasiparticle DOS 
instead of the real DOS. From (PJ it is easy to derive 
k(A)/k(0) = 1/(1 — A/c(0)£), in which all the correlation 
effects are contained in the purely electronic compress- 
ibility calculated at A = 0. This implies that, if the 



FIG. 2: The phase diagram for U/t = 10 as a function of 
doping S (n = 1 — 8) and A. 

Coming back to the results shown in Fig. 1, the n-/i 
curves of the metallic solutions have a small finite slope 
when they approach n = 1. However, if we consider 
the compressibility at fixed density as a function of A, 
it turns out that the e-ph interaction is substantially in- 
creasing the compressibility. In Fig. [31 we present the 
inverse of the normalized compressibility as a function of 
A at U/t — 5 for different values of the electron density, 
ranging from n — 0.70 to n — 0.95. An enhancement 
of k(X)/k(0) that varies from almost 6 in the n = 0.70 
case, to around 3 in the n = 0.95 case corresponds to the 
decreasing linear behavior of all the curves seen in the 
figure. We also calculated the values of the renormal- 
ization factor Z which is inversely proportional to the 
quasiparticle density of states. Since we found that Z 
varies much less than k as a function of A, such an en- 
hancement of the compressibility cannot be due to mass 
renormalization alone. This also indicates that, within 
the Landau Fermi liquid approach, the total Fq is ex- 
pected to be renormalized by the e-ph interaction like in 
the large— N calculations. We extract informations on 
this renormalization of the Landau amplitudes by simply 
assuming that a relation similar to Eq. holds also for 
the DMFT solution. We therefore fitted the DMFT data 
with the relation 
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*(A) = 

k(0) l-/5Art(0)r 



(4) 



where (3 is the only fitting parameter. The results of 
this fit, together with the values of k(A) in the inset, 
are shown in Fig. For almost the whole range of 
densities the fit is extremely accurate up to A = 2.2, with 
just a small correction to the ideal large-N result (3=1. 
The closer we get to half-filling the more f3 approaches 1, 
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solution with respect to the metal. This effect, besides 
a quantitative modification of the coexistence region be- 
tween the two solutions, determines a first-order transi- 
tion between an insulator and a metal at some value of 
the chemical potential fj, c , and to a PS between the Mott 
insulator with n — 1 and a metal with a finite density 
n ^ 1. This PS occurs for A ~ 1. For substantially larger 
values of the e-ph coupling a polaron crossover occurs 
very close to half-filling and PS takes place between this 
polaronic state and a non-polaronic metal. 

We acknowledge fruitful discussions with S. Ciuchi and 
support from MIUR Cofin 2001, prot. 2001023848, and 
from INFM, also through PA_G0_4 project. 
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FIG. 3: Normalized inverse compressibility [k(A)/k(0)] _1 as a 
function of A at U/t = 5 and for different densities, compared 
with the fit described in the text. Inset: Compressibility vs 
filling for various values of A = 0.0,0.8,1.2,1.6,1.8,2.0,2.2 
(from the bottom to the top). 



even if at n — 0.95 the strong-coupling data display some 
fluctuations. The values of (3 are reported in the figure. 

In light of the ability of the simple Eq. (0J to de- 
scribe the way the e-ph interaction modifies the scenario 
determined by electron correlation, the small values of 
A needed to obtain a divergent compressibility in the 
large— N approach 0, E3| can be related to the smaller 
value of the electronic Fq found in the large— AT with 
respect to the one found in the present DMFT approach. 

Our findings could be related to some experimental 
observations. First of all the presence of a (phonon- 
induced) PS close to a Mott insulating phase is suggestive 
of a chemical potential plateaux inferred from photoemis- 
sion experiments in underdoped cuprates 21] . The for- 
mation of doping-induced metallic states with a chemical 
potential slowly varying with doping inside a Mott gap 
was also recently deduced from tunneling experiments 
|22| . In the same experiment, the small mixing of the 
metallic (superconducting) phase with the antiferromag- 
netic insulating one indicates that a spatial separation of 
the two phases can have small interfaces suggesting that 
the presence of long-range Coulomb forces can easily in- 
duce small-scale (presumably stripe-like) domain forma- 
tion. We also remarkably found that for rather large 
values of the e-ph coupling (A > 2 in the U = lOt case) a 
PS can occur between an insulator away from half-filling 
with a finite density of polarons and a good metal with 
nearly free carriers. Also this situation of coexisting (but 
spatially separated) polarons and free carriers could find 
a realization in underdoped cuprates |23| . 

In conclusion, we have studied the effect of the e-ph in- 
teraction on the density-driven Mott transition for large 
values of U. The e-ph interaction favors the insulating 
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